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We compute, for the first time using lattice QCD methods, radiative transition rates involv- 
ing excited charmonium states, states of high spin and exotics. Utilizing a large basis of in- 
terpolating fields we are able to project out various excited state contributions to three-point 
correlators computed on quenched anisotropic lattices. In the first lattice QCD calculation of 
£NJ ■ the exotic 1 h ?y c i radiative decay, we find a large partial width r(?y c i — > J/tpj) ~ 100 keV. 

[ We find clear signals for electric dipole and magnetic quadrupole transition form factors in 

• Xc2 — * J/ipJ, calculated for the first time in this framework, and study transitions involving 

Ph ! excited tp and Xci,2 states. We calculate hindered magnetic dipole transition widths with- 

£f-^ • out the sensitivity to assumptions made in model studies and find statistically significant 

signals, including a non-exotic vector hybrid candidate Yhyb? — ► VcJ- As well as comparison 
to experimental data, we discuss in some detail the phenomenology suggested by our results 
' f-H ' | and the extent to which it mirrors that of quark potential models and make suggestions for 

O • the interpretation of our results involving exotic quantum numbered states. 
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I. INTRODUCTION 

£S| . The charmonium system is often described as the "hydrogen atom" of meson spectroscopy. 

Being reasonably non-relativistic, it is explained fairly well by potential models, at least 
CN ' below the open-charm (DD) threshold. More recently it has also been studied using effective 

field theory approaches (such as pNRQCD) and QCD sum rules. There has lately been a 
resurgence of interest in the charmonium system with the i?-factories, CLEO-c and BES 
finding missing states, making more accurate measurements of properties of these states and 
discovering a number of new resonances that are not easily explained by the quark model. 
This has spurred renewed theoretical interest with much speculation as to whether these 
states are hybrids or multiquark/molecular mesons. To date there are no charmonium states 
having manifestly exotic J PC such as 1 ^,0^ , 2 H , that would directly signal physics not 
present in potential models. 

The states below open-charm (DD) threshold can not decay via an OZI allowed strong 
decay and so have reasonably narrow widths. Their radiative transitions can therefore have 
significant branching ratios and are experimentally accessible. The transitions from and 
production of the, as yet unobserved, exotic 1 ^ are particularly interesting. A lattice cal- 
culation of transition form factors of excited charmonia is therefore timely and this is the 
first such study. The corresponding excited charmonium spectrum was calculated in lattice 
QCD in Ref. Transition form factors of the lightest few charmonia, those ground states 
accessible with interpolating fields f/T-i/;, were calculated in Ref. 0]; this work brought to the 
attention of CLEO-c experimentalists the discrepancy between the lattice calculated value of 
T(J/ip — > r/c-f) (and indeed the values predicted in most model calculations) and the single 
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experimental measurement of this from Crystal Ball 0]. In a tour-de- force analysis 0], a much 
more reliable value was extracted from CLEO-c data that is in much better agreement with 
theoretical estimates. 

The calculation we will present is performed in the quenched approximation, neglecting 
altogether the effect of light quark degrees of freedom. As such it is rather directly related 
to the simplest quark potential models in which charm quarks move in a static potential of 
assumed gluonic origin. Attempts have been made to add in the effects of light-quark loops 
to these models^, [a, 0], in some cases finding that these effects can be largeQ|. We will 
address this possibility in light of our results. 

A strong motivation for developing the lattice QCD techniques required to extract excited 
and exotic state radiative transition matrix elements is the versatility of the method. We can 
use these methods, tested here in charmonium, at any computationally feasible quark mass. 
In particular this opens up the possibility of computing in a framework close to QCD the 
meson photocouplings that appear in the meson photo-production process to be utilized in 
the JLab 12 GeV GlueX experiment 0]. In this paper we will also compare results with the 
flux-tube model of gluonic excitations which to date is the only theoretical guide to the size 



of the hybrid couplings and hence production rates[lQ. 111!]. 

The paper is structured as follows: We begin in section [II] with a description of the 
technology used to construct three-point correlators and to project on to the contribution 
due to various excited states. In section IIIII we present our results for the transition form- 
factors between various meson states. A discussion of the phenomenology of these results in 
terms of experiment and models follows in section IIVI before we conclude in section |V] Two 
appendices which consider some technical details complete the manuscript. 



II. TECHNOLOGY 

In this paper we will explore radiative transitions with charm-mass quarks using the 
quenched anisotropic lattices and clover fermion action described in [l|]. Radiative transition 
matrix-elements follow from vector current three-point functions, whose construction using 
sequential source technology is described in 1- As in i, only connected diagrams are 
considered with the assumption that disconnected diagrams are negligible in comparison. 
The new item of technology in this study is the use of a large basis of meson interpolating 
fields as explored in two-point functions in [l[ with further interpretation in [l2l |. 

These operators are lattice-discretised versions of the gauge-invariant fermion bilinears, 

i;(x)T*D l *D j ...^(x), 

where D = D — D is the gauge-covariant derivative operator. In order to improve the overlap 

on to lower lying mesons, the quark fields may also be smeared over space as e4 CT DiDi t(;(x). 
Linear combinations of these basic operators have been constructed which transform irre- 
ducibly under the lattice cubic symmetry and which have rather simple interpretations in the 
limit of zero lattice spacing, for details see [H, Q3] ■ The operator set is large enough to have 
a considerable redundancy within a given quantum number sector, which can be utilized to 
extract excited states. 
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A. "Ideal" operators and eigenvector projection 

In [H the spectrum of charmonium was extracted by solution of the generalised eigenvalue 
problem 1 

Cij(t)vJ = A n (t)Cy(ioK (1) 
^ nt C ii (ioX = <5 nm 

which constitutes a best solution for the spectrum in the variational sense. The matrix of 
two-point correlators is constructed as Cy(i) = (0| Y^g 0) 1 0) where the sum over 

lattice sites forces the three- momentum of single particle states to be zero 2 . This method 
relies upon a redundancy of operators Oi in any given quantum number sector. 

Fitting the time-dependence of the principal eigenvalues, A n (t), gives us the masses (or 
energies at finite three-momentum) of the states with the quantum numbers of the hermitian 
operators Oi t j. Details of the two-point analysis, including the importance of the choice of 
timeslice to, are given in [if. The interpretation of the eigenvector v n in the solution of Eqn. Q] 
can be expressed as follows: weighting the operators by this vector gives the optimal operator 
within the limited operator space for the state labeled by n. It is convenient to normalise these 
"ideal" operators in a manner which accounts for the value of t : 0" = ^j2E~ n e- EnU, l 2 v\Oi. 

Note that these eigenvectors are trivially related to the vacuum-operator-state matrix 
elements or "overlaps", Zf = (n|0j|O), that appear in a bound-state spectral decomposition, 

yn* yn 



by Z?vf = y/2E~ n e E " ta l 2 5 ntm or 



Z^=y/2E^^' 2 vfC ij {t Q ). (2) 

The procedure we shall follow in this paper is to compute three-point correlators hav- 
ing at the source (located on a fixed timeslice, tj) a (smeared) local operator of the form 
?/>(0, ti)Tip(0, ti) - all possible gamma matrices can be considered for the computing cost of 
a single "forward" propagator. In this work we consider V = 75, 7«, 1 giving access to pseu- 
doscalar, vector and scalar states at the source. 

At the sink (located on a fixed timeslice, tf) we use sequential-source technology using a 
broad selection of local and derivative-based operators as described in [l| - for each operator 
at a given momentum (usually pf = 000) we have the computing cost of a single "backward" 
propagator. At the vector current insertion (inserted on all timeslices tj < t < tf) we insert 
the local vector current with all possible lattice three- momenta q up to \q\ 2 = 4. Translational 
invariance ensures momentum conservation, selecting the correct value of pi out of the sum 
over all momenta produced by a local operator. 

Cth&w kAtf) = (o| e- ipr %&t f ) ■ *) ■ M U) 0) (3) 



1 repeated indices are summed 

2 In some cases we also considered non-zero three-momentum correlators - this will be discussed later 
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Three-point correlators for "ideal" operators are constructed from the "raw" correlators, 
eqn[3]by projecting with the appropriate eigenvector 3 , 

C r ^{pi,Pf,U,t,t f ) = (o\j2e~ ifr ^ n &tf) ■ ^2^%(y,t) ■ $@,U)n/>&U)\o 

z y 

= ^/W n e- E ^l 2 vfC riij {p u p } -t h t,t f ). 

The net effect is to provide us with a three-point correlator that, while it has multiple 
states contributing through the local source operator V>(0, ti)Tip(0, tj), should have only a 
single state, n, contributing from the sink operator Q n . In principle we could improve this 
further by using the full basis of operators at the source as well as the sink, but using 
sequential source technology this quickly becomes rather expensive 4 . 

In fact we can be a little more precise about the statement that only a single state con- 
tributes at the sink - in considering the analysis of the two-point correlators we determined 
that to can be considered to be the timeslice on which the correlator matrix (of dimension 
N) is saturated by N states - closer to the operator we require more states to describe the 
correlator matrix [l|, Q3] ■ It follows that if we are to timeslices away from the sink operator 
position tf we can be quite confident that we have contribution only from the single state 
labeled by n. When we come to fit the time-dependence of the three-point correlator we will 
not use timeslices any closer to the sink than t = tf — to- 

In figured] we show the effect of the eigenvector projection on to "ideal" states in the case 
of a pseudoscalar source operator ^7 ip and a selection of pseudoscalar sink operators. We 
clearly see that the projection on to the ground state "ideal" operator produces a considerable 
"flattening" of the correlator moving toward the sink. We do see some curvature beginning 
within to = 8 timeslices away from the sink as expected 5 . More importantly for this study we 
see that projection on to the "ideal" operator for the first excited state yields a clear non-zero 
signal from which one can extract a transition matrix element. 

There is in fact a slight technical subtlety to be dealt with in the eigenvector projection on 
to "ideal" operators that arises from the method of solution of the two-point problem. Solving 
the eigenvalue problem gives eigenvectors on each timeslice, v*(t2 P t) - if the solution is to be a 
true spectral representation there should be no time dependence in the eigenvectors, at least 
for t2pt > *o where the correlator matrix is supposed to be saturated. In |l| we saw that over 
a considerable range of t2 P t, the Z values (obtained trivially from the v, see equation [2]) were 
flat. We explore the effect of any such possible time-dependence on the three-point functions 
by performing the projection for eigenvectors belonging to a range of timeslices, v^(t2 P t)- In 
figure [2] we show for a ground state and first excited state projection the *2pt dependence 
on the projected correlator. It is clear that there is no strong dependence on t2 P t and we 
choose to average over the projections for various t2 P t to reduce configuration-by-configuration 
fluctuations that may be caused by the generalized eigensystem solver. 



3 A somewhat different appiication of the same basic idea has recently been presented in 
In addition analysis of these correlators would require detailed knowledge of the finite-momentum behaviour 
of the derivative-based operators - we will consider this in appendix A and find that it is not entirely trivial 

5 Indeed if we fit the time-dependence of the curvature at the sink with a single exponential we find that it 
corresponds to a mass heavier than the heaviest state extracted from the eight operator two-point correlation 
matrix 
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FIG. 1: Blue and yellow data - suitably normalised pscudoscalar- vector current-pscudoscalar cor- 
relators for a set of different sink operators. Red data - correlator formed by projection with the 
ground state eigenvector vf\ Green data - correlator formed by projection with the first excited 
state eigenvector v!> 




FIG. 2: Eigenvector projected correlators as a function of current insertion time t. For each t value we 
show the result of projecting with each of 1 (t^pt) for 11 < t2 P t < 22. Also shown is the average 
over the various t2 P t values. 



I. 



B. Transition form- factors 

The "ideal" operator-projected three-point functions are related to the transition matrix 
elements between states as follows: 

Cvfin{pi,Pf,ti,t,tf) = 

. \ p- E n{tf~t) . . p -E m (t-ti) , . 

n»(0,0) hipf)) ^ ^(/n(p» ^(0,0) i m (pi}) ^ (»m(Pi) ^(0,0) O) (4) 



2£, 



m 



While all possible eigenstates with the quantum numbers of Oi will contribute at the 
source (the sum over m), the eigenvector projection ensures that we need only consider 
the single state n at the sink. Our normalisation of O n ensures that if / is spin-zero, 
(o|fi n (d,0)|/ n (p>)) = 2E n . 

The Minkowski-space transition matrix element (/n(p/)|j M (0, 0)\i m (pi)} can be decom- 
posed in terms of multipole form-factors multiplied by Lorentz covariant combinations of 
the momenta and (if appropriate) polarization tensors of the particles labeled by i, f: 
^2 k Fk(Q 2 )K^(pi,pf,ei,e*j). The general technique for obtaining these covariant multipole 
decompositions is given in the appendix of [3]. Inserting this decomposition into eqn. IH 
and performing the implicit sums over helicity (in the case of spin > 1), the three-point 
correlator on each timeslice, t, can be expressed as a linear sum of known "kinematic" and 
"propagation" factors times the unknown multipole form- factors 6 , 

,m,n{Q )• 

m k 

In practice we opt to solve this linear system 7 on each timeslice, t, for a set of "effective", 
t-dependent form-factors, i 7 fcn(Q 2 ,t), using the propagation and kinematic factors for the 
ground state at the source(m = 0): 

Crtin(Pi,Pf, k,t, t f ) = P , n {Pi,Pf,ti, t, t f ) ^2 K kfi,xSP^Pf) F kAQ 2 

k 

The propagation factor Po in (Pi,Pf, U,t,tf) contains the overlap Z^(pi) and the energy E${pi). 
These are determined in the analysis of two-point function at finite momentum. Since the 
state masses at zero momentum are rather precisely determined by the variational solution, 
we use the continuum dispersion relation to obtain E(p) = \Jm? + \p\ 2 . Z values then 
follow from a linear fit to the finite-momentum two-point correlators supplying the known 
exp (—E n (p)t) factors. 

The extracted "effective" form factors have the property that Fk jn (Q 2 ,t) « F^ ^Q 2 ) + 
O (e~( £ ' m = 1- - Bm =°)(* - * 4 )) so that away from the source the excited state contributions die away 
and the form-factor plateaus to the i m =o - > fn value. We fit the time-dependence of these 
extracted "effective" form-factors to a sum of exponentials where the energy dependence is 
that extracted from spectrum studies 8 . A typical example is shown in figure [3J We retain 



6 This is essentially the same decomposition presented in Q| 

7 by jackknifed SVD as described in 

8 With the exponential dependence supplied this is a linear fit for the coefficients 
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FIG. 3: "Effective" time-dependent form- factor. Fit with constant plus a single exponential shown in 
red, constant fit value shown in green. Fit is linear for the coefficients since the excited state energies 
appearing in the exponent are determined from solution of the two-point function problem. 



only the plateau value, discarding the relatively unreliable source excited state information. 

The expansion of transition matrix elements as performed above relies upon the continuum 
Lorentz symmetry which is broken on our cubic lattice. A rigorous study would consider the 
expansion in terms of irreducible representations of the cubic group at finite momentum 
and would typically involve a larger number of form-factors, some of which would tend to 
the multipole form- factors as a — > and some of which would have to vanish. We have 
not performed such a decomposition but suspect that it is not important at this lattice 
spacing since when we solve the correlator-form-factor linear system assuming continuum- 
like decompositions we typically obtain \ 2 P er d.o.f. very close to 1, suggesting that there is 
little need to enlarge the basis space. 

The fact that charmonium states are eigenstates of charge conjugation (C), coupled with 
the photon having C = — 1, means that all charmonium radiation transitions are between 
states of opposite C. The origin of this is that the photon couples equally to the charm quark 
and the charm anti-quark. In our lattice calculation we need not impose this symmetry; by 
coupling only to the quark we can obtain transitions between states of equal C. One way 
of viewing this is that it is like having a u, d pair but as heavy as the charm quark - this 
can useful for comparison with models and we will show results for a number of C-violating 
transitions. 



Strictly speaking the kinematic factors for excited m states need not be proportional to the ground state 
values such that excited state multipoles may "leak" into the wrong ground state multipole - clearly this 
effect will fall off for t^ti. An alternative approach is to include the excited state terms directly into the 
decomposition fit, expanding the space to also be over timeslices, the time-dependence being the discrimina- 
tor for the various states contributing at the source - this guarantees the right kinematic factors for excited 
states. In all cases to be presented the analysis was cross-checked using this fitting method (the excited 
state energies obtained from applying the dispersion relation to the excited state masses reliably extracted 
using the variational method) and agreement within statistical fluctuations found. 
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level 


mass / MeV 


suggested state 


model assignment 





3106(2) 


J/il> 


l 3 5i 


1 


3746(18) 


V>'(3686) 


2 3 ^ 


2 


3846(12) 




lat. artifact 


3 


3864(19) 


V>"(3770) 


1 3 D 1 


4 


4283(77) 


^('4040') 


3 3 Si 


5 


4400(60) 


Y ? 


hybrid 



TABLE I: T l spectrum extracted from two-point functions, suggested experimental state analogues 
and quark model bound-state assignments made in [12|. 



III. RESULTS 
A. Vector current renormalization 

Since we use the local vector current in the construction of our three-point functions we 
need to determine the renormalization constant, Zy, to relate the extracted matrix elements 
to physical matrix elements. We do this by insisting that the pseudoscalar form- factor at 
zero Q 2 takes the value 1. On an anisotropic lattice we should allow there to be differ- 
ent renormalization constants for temporally and spatially directed currents, indeed we find 
Zyf s ) = 1-23(2) and Z v ^ = 1.118(6). In all results presented below only the spatially di- 
rected current is used. See appendix B for a discussion of the effect of improvement of the 
vector current. 



B. Scalar - Vector transitions 

The first results we will present concern transitions between scalar (0 ++ ) and vector 
(1 ) states of charmonium. Using a quark-smeared operator ■ipip at the source and 13 
vector operators 10 at the sink we extracted transitions between the ground state scalar XcO 



and the lowest six vector states. In table U we show the spectrum reported in [12|], possible 
comparable experimental states and a model-dependent state assignment. 

This transition is characterized by two multipole amplitudes, a transverse electric dipole 
E\(Q 2 ) and a longitudinal C\(Q 2 ), the matrix element decomposition being 



(S(p S )\f(0)\V(py, A)) = fT 1 ^ 2 ) E l {Q 2 ) \n(Q 2 )e^p v , A) - e(p v , A) • p s {p v Pv ■ Ps - m 2 v p^ s ) 



+ Cl ^%' ^ mve(p v , A) - ps Pv PsiPv + PsY ~ m 2 s p v - m v p^ 



(5) 



where £l(Q 2 ) = (py ■ ps) 2 — m 2 s my. 

The "vector" operators used are, in fact, all in a particular irreducible representation of 
the cubic group at rest, namely T\ which in the continuum limit contains spins 1,3,4.... By 



In the notation of [l[ they are the quark-smeared versions of 7i,7i7o,ao x Vti,<ii x Vti,P X Dti,P(2) x 
Dti,7t x Bti,7T(2) x Bti and the unsmeared versions of 7^,7^70,01 X Vti,tt x Byi,7T(2) x Bti. In the 
Pf = (100) case the p( 2 ) x D and 7T( 2 ) x B operators were not included. 
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sink suggested 
level transition 


W A/GeV" 2 


r lat /keV r oxpt /keV 


x c o -> J 1^1 

1 if — > Xc07 

3 V" -> Xco7 

5 Ihyb. -> Xc07 


0.127(2) 409 < 12 ) 
V ; 1.14(5) 

092(191 164(55) 
[fixed] 

0.265(33) 324 f) 
v ; 0.58(56) 

0.00(3) linGar 
fit 


199(6) 131(14) 

265(66) 199(26) 
< 20 



TABLE II: Results of fit to lattice data using equation [5] Partial decay width computed using 
fitted value of -Ei(O) and physical phase space (where known). All errors are just lattice statistical. 
Experimental partial decay widths from [3, [l6| . 



considering the degeneracy pattern and the values of the overlaps Z • , the two-point function 
analysis strongly suggests that the 2nd excited state in the T x sector is in fact a spin-3 state 
and as such we do not present results involving this state here [l|, 03]. Additionally we do not 
show the amplitudes for the 4th excited state 11 as the signals are statistically consistent with 
zero for all Q 2 . Shown in figure 0] are the electric dipole transition form- factors for the four 
relevant states. See appendix A for a discussion of the inclusion of pf = (100) correlators in 
this analysis. 

In each case the Q 2 dependence of the lattice data has been fitted with the form, 

_ Q 2 

Ei(Q 2 ) = £i(0) (1 + XQ 2 ) e (6) 

whose motivation is described in [2j. The results of the fits are given in Table M One point 
to make here is that in our lattice calculations, the photon only couples to the quark and not 
the antiquark and we do not explicitly include the electric charge of the quark. Therefore we 
actually compute and present Fk(Q 2 ) which are related to the physical multipole amplitudes 
by Ff t (Q 2 ) = 2 x |e x F^{Q 2 ). Note also that the multipole amplitudes have mass dimension 
1 and hence we plot them in temporal lattice units, where a^ 1 = 6.05 GeV is determined 
from the static potential. Partial decay widths follow by averaging over initial helicities and 
summing over final helicities which gives, for a general A — ► B^f decay: 

2 J a + 1 9 m A t-f I 

Where we have experimental masses for states we use these to c omp ute the phase space, 
otherwise we use the value extracted from the lattice calculation [3. fl2jj. 

The ground state transition XcO ~^ J 1^1 form-factor shows behavior rather similar to 
that found in [3] which used the same lattices but Domain Wall fermions rather than Clover 
and which did not make use of multiple sink-operators and operator projection. We will 
consider this in a little more detail in Appendix B where we will consider the effect of 0(ma) 
improvement of the local vector current. We will discuss the results in comparison with 
experiment and with quark potential models in section HV Bl 



In [13] this state was identified with the quark model state 3 3 Si 



10 




0.5 1 1.5 2 


2.5 3 3.5 

QVGeV 








: 1 



\ 

□ 

A O 

, , , 1 , , , , 1 , 


p-000 

P r = i 'oo 
, , 1 , 



0.5 1 1.5 2 2.5 3 3.5 



QVGeV 



r^x co r 


□ 


p=000 




o 


p~100 


^.CLEO'06 












. ! . . . . i . . . . i . . . . i . . . . i . . . . i . . . . i . 


, , , 1 , 



0.5 1 1.5 2 2.5 3 3.5 

QVGeV 



hyh? ~ * %c0 y 




QVGeV 



FIG. 4: Electric dipole transition form-factors Xco i>- Plotted is the form-factor in temporal lattice 
units against the photon virtuality in GeV 2 . Fits to the lattice Q 2 dependence as described in the text. 
Experimental points at Q 2 = are extracted from experimental decay widths taken from 



11 



sink 
level 


suggested 
transition 


V{U) 


/3/MeV 
A/GeV- 2 


r lat /keV r oxpt /keV 


U 


j t n i. 
J IV 




l.oy^o ) 


513(7) 
[fixed] 


2.51(5) 


l.o5(29 j 




v'- 


%7 




530(110) 
4(6) 


0.4(8) 


0.95(16) 
1.37(20) 


3 




-> »7c7 


0.27(15) 


367(55) 
-1.25(30) 


10(11) 




■5 


^hyb. 


->• ^7 


0.28(6) 


250(200) 
[fixed] 


42(18) 





TABLE III: Results of fit to lattice data using equation [|5J Partial decay width computed using 
fitted value of V(0) and physical phase space (where known). All errors are purely lattice statistical. 
Experimental partial decay widths from [J, [l5| . 



C. Vector - Pseudoscalar transitions 

Using the same T\ operator set at the sink and the quark-smeared ip^ 5 ip operator 
at the source we obtained results for the single magnetic dipole form-factor in the vector- 
pseudoscalar transition. The decomposition used 12 is 

(P(p P )\f(0)\V(p v ,X)) = W{ f ] e^pp Q p Vf3 e,(p v ,X). 

rap + my 

Figure [5] shows the form-factors for the transition between the lightest four vector states 
(ignoring the suspected 3 intruder and the noisy ■0('4O4O') state) and the r/ c . The fit-form 
in eqn [6] was again used, with the fit results presented in table IIIII We refer the reader to 
the paper [3] for a discussion of the systematic error introduced into the phase space by an 
inaccurate lattice estimate of the hyperfine splitting. 

D. Axial - Vector transitions 

In transitions between axial (1 ++ ) and vector (1 ) states there are two transverse form- 
factors, electric dipole (Ei) and magnetic quadrupole (M2), and one longitudinal form- factor 
(Ci). The matrix element decomposition takes the form 

(A(p A ,X A )\f(0)\V(p v ,Xv)) = ^hm^ upa (PA-Pv), 

Ei(Q 2 )(p A +p v )p(2ra A [e*{p~ A ,X A ) • Pv]^u(pv, Ay) + 2m v [e(pvAv) ■ PaK(pa, Xa) 

+ M 2 (Q 2 )(pa + Pv) P (2rn A [e* (p~a, Xa) • Pv]^(pv, Ay) ~ 2my[e(py,Ay) ■ p A ]e* u {p A ,X A ) 
An(Q 2 )el(p A ,X A )e p (p v ,X v ) 



+ {PA+Pv) P (m A - m\ + q 2 )[e*(p A ,X A ) ■ p v ] e v {p v ,X v ) + (m\ - m 2 v - q 2 )[e(p v , X v ).p A ] e*(pl, \a) 



note that this is not the conventionally normalized magnetic dipole amplitude 
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FIG. 5: Magnetic dipole transition form-factors ip — ► t] c . Plotted is the dimensionless form-factor 
against the photon virtuality in GeV 2 . Fits to the lattice Q 2 dependence as described in the text. 
Experimental points at Q 2 = are extracted from experimental decay widths taken from 0, . 
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FIG. 6: (a) Electric dipole and magnetic quadrupolc form-factors for the transition Xd ^ J/^l- 
Q 2 dependence fitted with equation [5] and extrapolated to the physical photon point Q 2 = for 
comparison with experimental data from [lij]. Relative sign of E\ to M2 is relevant, (b) Ratio of 
magnetic quadrupole to electric dipole form-factors. Colored bands represent fits with various fit 
functions. 



In this case we computed with the vector at the source, produced using quark-smeared 
•tjyyiip, and a set of ten axial (actually T^~ + ) operators at the sink 14 . In figure [6ja) we 
present the two transverse form-factors for the ground state axial meson and in figure [6](b) 
their ratio. The amplitudes were each fitted with a form like eqn. [6] which unfortunately 
is not as well constrained as in the scalar case owing to the kinematic factors preventing 
a slightly time-like Q 2 point corresponding to Pi = p) = (000). The results of the fits are 
atEi(O) = 0.23(3), p E i = 440(40) MeV, X E i = 0.71(30) GeV" 2 and a t M 2 (0) = -0.020(17), 
Pm2 = 450(50) MeV, Xm2 = 5(6) GeV -2 . This corresponds to a partial decay width of 
r(Xei - ► J/ip'y) — 270(70) keV which is in reasonable agreement with the PDG's average of 
320(25) keV. 

The ratio j^-(Q 2 ), shown in figure [6^b) , was fitted with various functional forms shown 



We could have used at the source to produce the axial meson, this operator at finite momentum 

also has overlap with the pseudoscalar state which severely limits its usefulness. 

7i75, p x Vti, P(2) x Vti, a,i x Dti, bi x Bti in both quark-smeared and unsmeared versions. 
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FIG. 7: Electric dipole and magnetic quadrupolc form-factors for the transition x'd ~ * J/i^l- Q 2 
dependence fitted with various forms shown by the colored bands. Relative sign of E\ to M 2 is 
relevant. 



by the shaded bands yielding ^(0) = —0.20(6) where the error includes a crude estimate 
of the systematic error due to the uncertainty in fitting form. The ratio of the extrapolated 
values from the separate fits to Ei, Mi gives ^(0) = 023(3)'^ = — 0.09(7). Clearly without 
data points at smaller Q 2 or some certainty about the expected Q 2 dependence, we cannot 
constrain this any further and hence cannot make a particularly meaningful comparison with 
the PDG average — MM _ = -0.002+R H?f . 

& y/Ei(0) 2 +M 2 (0) 2 -°- 017 

Form-factors for the transition from the first excited axial state, x' c i down to the J/ip 
are shown in figure [7] where multiple fit forms were used, all returning a x /-^dof close to 
one. The estimates for the physical photon point thus obtained are at-Ei(O) = 0.050(15) 
and <XfM2(0) = —0.004(14) where again we include a crude systematic error estimate for the 
fit-form variation. The E\ transition corresponds, for a x'd & t 4.1 GeV, to a partial decay 
width r(x' cl -» J/ilry) = 21(12) keV. 



E. Tensor - Vector transition 



The transition between the lightest 2 ++ state and the lightest vector state was not con- 
sidered in [3] as that study used only local fermion bilinears ('ijj(x)F'ip(x)) to produce states 
- a spin-2 particle cannot be produced by any such operator. Here we use a set of six opera- 
tors projected into T^~ + and E ++ irreps at the sink 15 . The multipole decomposition for this 
transition takes the following form, 



p x V, p( 2 ) x V, b\ x B in both quark smeared and unsmeared versions. 
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In order to have a constrained linear system for extraction of the five form-factors we are 
required to consider all five spin-2 helicities, which with respect to the lattice cubic symmetry 
are distributed in T2(3) and F(2) irreps. Strictly speaking these are independent irreps on the 
lattice and we should be careful about combining them, but any lattice symmetry breaking 
must vanish as a — > and in [l| we found that there were strong signals that the continuum 
rotational symmetry was restored to a good approximation already at this value of a. For 
example, for the ground states in the T 2 , E channels we have a high degree of mass degeneracy, 
and the values of the overlaps (Zt 2 /Ze) are compatible at the 1% level. Details of this analysis 
are in [lj]. 

The extracted form- factors and fits to the Q 2 dependence are shown in figure [8] where the 
expected hierarchy |Fi(0)| > |M2(0)| 3> |F3(0)|, to be discussed in section HVB| is observed. 
The description of the lattice three-point correlators by the continuum decomposition, eqn 
[7] (which does not take any account of T2/E discretization differences) is excellent, with 
typically x 2 /-^dof ~ 1- As a simple systematic test of the degree to which T2/F discretization 
differences could affect the determination of the form-factors we tried deliberately introducing 
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FIG. 8: Electric dipolc, magnetic quadrupole and electric octopole form-factors for the transition 
Xc2 — > J/^Pj- Q 2 dependencies fitted with equation El an d extrapolated to the physical photon point 
Q 2 = for comparison with experimental data from |15| . Relative signs are relevant. 



such a difference by arbitrarily multiplying all E correlators by a factor. Even with a 10% 
increase in the magnitude of all E correlators, the form-factor values changed by less than 
the statistical error, while the x 2 /-Wdof increased to around 6. If we where to see such large 
X 2 in our raw data (which we never do) we would not trust the results - we conclude that our 
form- factor extractions are not being strongly affected by TijE discretization differences. 

The amplitudes at Q 2 = correspond to a partial decay width of T(xc2 — ► J/tp"f) = 
380(50) keV to be compared with the PDG average of 406(31) keV. The electric dipole fit 
parameters are P = 550(80) MeV and A = -0.39(1) GeV" 2 . The fits yield values for the ratios 
of multipole amplitudes of m 2 (0) _ _q ^q^\ j e 3 (o) _ 

L ^ 7^1(0)2+^/2 (0) 2 +£3(0)2 V ' (O) 2 +A/ 2 (0)2+^3 (0) 2 

0.010(11) which in the case of the first ratio is considerably larger than the PDG average of 
—0.13(5) but which does appear to be of the correct sign. 

We also extracted the x' C 2 ~^ J/^l transition form- factors 16 as presented in figure [9](a). 
The large value of £"3 might be surprising when compared to the ground state result - a 
number of simple systematic tests were performed to investigate if this might come about 
through various lattice effects with the result that we were unable to change the values outside 
statistical error bars by any reasonable adjustment. The amplitudes at Q 2 = correspond to 
a partial decay width of 20(13) keV. In figure [91(b) we present the transition form-factors for 
the next excited tensor state 17 , x% ~ * J/^l where the hierarchy |iTd(0)| 3> |M2(0)|, ^(O)) 
appears to be restored and where we predict a partial decay width of 88(13) keV. In section 
IIVBI we will propose a simple explanation for all these observations in the framework of a 
non-relativistic quark model. 



F. Exotic transitions 

Our principal focus here is the exotic 1 ^ state, r/ c i, here 18 found to be at 4300(50) MeV. 
On a quenched lattice, provided we can eliminate the possibility that the lightest state in 
is part of a non-exotic 4 ^ (see [l^j for support of this elimination), we can be fairly certain 



We found a mass of 4115(28) MeV for the x'c2 
We found a mass of 4165(30) MeV for the X c2 

But note our comments regarding the box-size as a cause of systematic error 
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that this state is a hybrid, having an excited gluonic field in addition to a charm-anticharm 
quark pair. This is in contrast to having a higher quark number Fock state which, since we lack 
light-quarks altogether in this calculation, could only arise for states having mass near 4m c 19 . 
A charge-conjugation allowed decay of this meson would be r) c \ — ► J/ip'y having transverse 
magnetic dipole and electric quadrupole multipole contributions. The decomposition of the 
vector current matrix element between two non-identical vector mesons takes the following 
form: 



and in a quenched calculation they would arise in a unitarity-violating way U\ 
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FIG. 10: Magnetic dipole and electric quadrupole form- factors for the transition T] c i(l + ) — > J/tpj- 
Point at Q 2 = 0.06 GeV translated to the physical photon point Q 2 = 0. 
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where IF = (p' + p)" - *2- 

We used a set of eight operators at the sink to produce the T^ - " 1 " state 20 . Our results 
projected on to the ground state in + are shown in figure [10] - the lattice data are not 
fitted as a function of Q 2 as very little extrapolation is required to associate the point at 
Q 2 = 0.06 GeV with the real photon point. The value of iWi(0) corresponds to a partial 
decay width T(rj c i — > J/tpj) = 115(16) keV. Note that this is no different in scale to many 
measured conventional charmonium transitions. 



In the nomenclature of [J] they are smeared and unsmeared versions of do(2) X Vti , &i x Vti , P x Bn , P(2) X 
Bti- The B operators induce an essential gluonic component through a factor proportional to F^ v in the 
continuum limit, the V operators reduce to covariant derivatives and hence have a factor A 11 . 
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FIG. 11: (a) — > Q~(+)~{ magnetic dipolc transition. Linear interpolation using the two lowest- 
lying Q 2 points gives the value shown at Q 2 = 0. (b) l - (+) — > + ' + ^7 electric dipole transition. Fit 
is with an exponential in Q 2 . 



If we also consider charge-conjugation symmetry violating decays (by coupling only to one 
quark) we have access — ► + ^ + ^7. These results can be compared 
to models of hybrid bound-state structure as we will do in the following section. Results are 
shown in figure [TTJ 

Another exotic J PC to which we have easy access is the 0" 1 state which appears in the 
A\ channel at 4465(65) MeV. We have computed a charge-conjugation symmetry violating 
decay of the lightest state with these quantum numbers to the conventional ground state 
vector 1~(~). This transition matrix element has the same decomposition as the vector- 
scalar in equation We used a set of three operators 21 to create the + ( _ ) state; in figure 
1121 we show the electric dipole form-factor for + (~) — ► 1 _ ( _ ^7. 



21 quark smeared i/j^tp and smeared and unsmeared ai x Bai 
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FIG. 12: Electric dipole form-factor for the transition — ► 1 — ('7. Fit is an exponential in Q 2 . 



G. Other assorted transitions 



We also computed some C-violating transitions which have allowed analogues in the light 
quark sector 22 , 1+C+) — > _ ^ + ^7 and 2 + ( + ) — > _ ^ + ^7. The first of these has a multipole 
decomposition identical to eqn[5]while the second takes the form 



(P(pvW(0)\T(pr, A)) = M 2 (Q 2 ) V2-^e^p Tp p Pa e UT (p T , \)p P 

Results are shown in figure [T3l A discussion in terms of non-relativistic quark models follows 
in the next section. 



IV. PHENOMENOLOGY 



The results presented in the previous section can be compared with experimental mea- 
surements and additionally with the predictions of various models describing heavy quark 
bound states. A rather successful approach to describing charmonium has come by consider- 
ing heavy charm quarks to be moving non-relativistically (or nearly so) in a static potential 
motivated by QCD. These quark potential models have been used to compute spectra, radia- 
tive transition rates, and, in many-body extensions of the theory, hadronic decay rates (e.g. 

In their simplest formulation, the gluonic field plays no dynamical role and one 
does not have exotic quantum numbered states having excited glue (hybrids). In order to 
consider such states one can construct models making specific assumptions about the nature 
of the gluonic field and its excitations, examples being the flux-tube modelfi^], "constituent" 



gluon models and models based upon a many-body treatment of QCD in a physical gauge [2 11]. 
In what follows we will address to what extent our lattice QCD results inform these models. 

Firstly we should discuss the degree to which the approximations we have made in our 
computation introduce systematic error to our results. 



namely a* — > jr 7 and 



7T ± 7. 
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FIG. 13: (a) Electric dipole form-factor for the transition — > '+'7. (b) Magnetic quadrupole 
form-factor for the transition 2 + W — > _ ' + ^7. Note that relative signs between states are not relevant. 

A. Lattice Systematics 

This calculation has been performed within the quenched truncation - as such one way 
to view it is as a calculation of a version of QCD having just one heavy flavor of quark 
with the neglect of heavy flavor quark loops being justified by their large mass. In this 
regard this calculation is rather directly comparable to quark potential models which in their 
simplest form also neglect the effect of dynamical light quarks. We set the lattice scale in our 
calculation using the Sommer parameter which is itself related to the static quark potential 
and the charmonium spectrum, so here too the comparison is fairly direct. In figure [14] we 
show the static potential extracted on this lattice along with the phenomenological forms 
used in a selection of potential models. 

Clearly then our calculation does not include the effect of charmonium states coupling to 
multi-hadron states containing light quarks, e.g. DD. There have been suggestions within 
extensions of quark-potential models that including such physics can have a considerable effect 
on the spectrum 0, 0] and radiative transition rates 0j of charmonium. We will attempt to 
address this in the discussion to follow. 

In terms of direct comparison with reality, in which there are three flavors of quark lighter 
than charm, another failing introduced by the quenched approximation is the incorrect run- 
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FIG. 14: Lattice static potential along with phenomenological potentials used in 0, EI, 22 | 



ning of the strong coupling owing to the beta function containing the wrong number of 
quark degrees of freedom. We expect this to show up in terms of inconsistent scale setting - 
something that has been observed previously f23| . 

Our calculation has been performed at only one value of the lattice spacing; as is well 
known one expects there to be systematic shifts in quantities owing to the lattice discretization 
that go away as one approaches the continuum. These are suggested to be particularly serious 
for the discretization of heavy quarks. At a single lattice spacing we cannot accurately 
estimate the size of such effects, although our use of the improved Clover action and an 
anisotropic lattice 23 should help a great deal in the reduction. We present in Appendix B 
a limited study of improving the vector current insertion in the manner described in [24]. 
In general the improvement effects are not large, suggesting that scaling to the continuum 
should not overwrite our results. 

In the spectrum calculations using the same action on the same lattices 0, 03, the pos- 
sibility that the spatial volume is too small to comfortably house highly excited states was 
raised. 

One way to reduce the systematic error introduced by extrapolating from finite Q 2 to Q 2 = 
would be to consider utilizing twisted boundary conditions on the propagator inversions 
to get momenta and hence Q 2 values rather close to Q 2 = - this has proven successful in 
form- factor computations [2o|| . 

Modulo the caveats that we have raised, we remain convinced that the results we present 
likely represent a faithful description of the pattern of physics of radiative transitions in 
charmonium, including the properties of gluonic hybrids. Most importantly we have clearly 
demonstrated that the technology of projecting three-point functions using "ideal" operator 
eigenvectors obtained in two-point function calculations works well in giving us access to ex- 
cited state transitions. All the possible sources of systematic error can be addressed in future 
calculations using a set of sufficiently large dynamical lattices of various lattice spacings. 



Some comparison of the Clover action with the Domain Wall action on these lattices is given in flj 
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B. Conventional state transitions and the quark-potential model 



Ref. 12| used information extracted from two-point functions to identify excited charmo- 
nium states found in a lattice calculation with the corresponding states expected in the quark 
model. The transition form-factors can also be compared to calculations in quark-potential 
models. Our quenched lattice calculations offer a rather direct comparison with these po- 
tential models: in both cases loops of mesons containing light quarks have been ignored. As 
discussed in Ref. @, [26|, the results of quark models are sensitive to approximations, such as 
the choice of frame, which are not an issue in lattice calculations - we will not discuss these 
further here. In the following, for simplicity, we use harmonic oscillator wavefunctions and 
consider the rest frame of the initial meson. 

Within the simplest non-relativistic model of a meson emitting a photon of momentum 
q, we find the following transition form factors for non-radially excited states undergoing a 
change of orbital angular momentum of one unit: 



Ei(\q\ 2 ) oc 
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where (5 is the harmonic oscillator wavefunction parameter. Here the r-values follow from 
the coupling of quark spin and orbital angular momentum as 
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Transforming from \q\ 2 to the invariant virtuality of the photon (Q 2 ) in the rest frame of the 
decaying meson this corresponds to 
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Note that this is of the general form given in eqn. 



In the same model, we have for transitions involving both a change in orbital angular 
momentum and a quark spin-flip, the expression 
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This applies to the charge-conjugation forbidden decays l + W( 3 Px) — ► ^ + -*( 1 5'o)7 and 
l+H^Pi) l"(-)( 3 5i)7. In addition, this same form is predicted for the M 2 form- 
factor in the transitions Xc1,2( 3 -Pl,2) - ► J/^^Si)^/ and the charge-conjugation forbidden 
2 +(+)(3p 2 ) _> o-(+)( 1 S"o)7 and l+^^Pi) -> l-(-)( 3 5i)7- In terms of Q 2 this gives 



Ei(Q 2 ) oc 

which we note is also of the form of eqn[6j 
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note that a typographical error in Q| is corrected here 
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Within this model we find that E 3 (Q 2 ) = for Xc2^P , i) — > J/ip( 3 S-\)"/. In fact this result 
follows from the general assumption of a single quark transition [271] and does not depend 
critically upon the details of the model. For a transition J — > J' (where J and J' are 
respectively the initial and final meson spins) the allowed k (for Ej. or depending on k 
and the parity change) are in general \ J — J'\ < k < J + J' . However, if only a single quark is 
involved in the transition there is a further restriction on the values allowed. If the interacting 
quark in the initial (final) meson has total angular momentum j (f) then \ j— j'\ < k < j+f- 
Hence in general for an n 3 P2 — ► n' 3 S\ transition the allowed k are k = 1,2,3. However in 
the single quark transition assumption we have j = 3/2 — ► f = 1/2 and k is restricted to 
k = 1,2. Note that it is possible to have E% ^ if the tensor meson is a 3 F2 (j' = 5/2) 
state, in which case k = 2, 3 and E\ = 0. It is also possible to have £3 ^ in transitions 
involving hybrid, multiquark or molecular mesons where there are other degrees-of- freedom 
able to carry some angular momentum. 



1. i\) — r] c transitions 

Our result for T(J/ip — > 7^7) is in reasonable agreement with the new experimental value 
from CLEO-c 0], the potential model calculations of 0, the pNRQCD result of ^ 



and the QCD sum-rule estimate of [30j]. Similarly our result for the "hindered" M\ transition 
T(ip' — > r/cj) is consistent with the experimental values [3, 15] and Eichten et. al. 28], but 
significantly smaller than the values in Ref. [l8| which depend sensitively on the overlap of 
orthogonal wavefunctions with a small recoil factor breaking the orthogonality. 

Li and ZhaoQ] study the possibility of hadronic loop contributions in J/tp and ip' radiative 
decays to r/ c and rj' c . They find the addition of loop contributions involving pairs of D and 
D* mesons can bring the model values from Ref. [3] in to line with experiment. Most 
significantly, they find large loop contributions, of order 10 keV, to the if)' — ► ^7 transition. 
A cancellation between the quark model contribution of order lOkeV (from Ref. [l^]) and 
these loop contributions results in a small transition width of ~ 1 keV. This appears to 
present a problem with respect to our result, which being quenched does not include loop 
contributions, since there is no room for a large loop contribution when our T{ip' — > rj c ^) 
is compared to the experimental value. This might indicate that the loop contributions of 
Ref. [8j are overestimated and we have identified one possible source for this in the coupling 
9D*D-y used. The partial decay width T(D*° —> D ^) ~ 800 keV assumed in Ref. 0] appears 
to be rather large for an Mi radiative transition and, with the branching ratio from the 
PDG[li], implies a D*° total width of 2.1 MeV which is at the upper limit allowed by the 



PDG. For comparison, Close and Swanson[3l[, using a non-relativistic quark model predict 
T(D*° — > D°7) = 32 keV. The corresponding decrease in the loop amplitudes would remove 
the problem in comparison with our result. In addition, there are considerable uncertainties 
arising from the estimated hadronic coupling constants in the loop contributions which are 
not discussed in Ref. 0|. These uncertainties may be particularly important because of the 
delicate nature of the cancellation. 

We have also been able to extract a signal for the tp" — > r/c'j which in the standard inter- 
pretation of the ip" would be a 1 3 Z?i — > 1 So transition. In a harmonic oscillator wavefunction 
basis the form-factor has the same leading |gj 3 behaviour as 2 3 Si — > 1 S'o7 indicating that it 
should suffer "hindered" suppression just like the ip' — > r\ c ^ relative to J/ip — ► 7^7. The lattice 
data for the amplitude appears to be in line with this, within a large statistical uncertainty. 
Discussion of the next excited vector state will follow in section HV CI 
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2. x-ip transitions 



The Xc.j J/^l partial widths we obtain are dominated by the electric dipole compo- 
nent and compare favorably with the experimental values. The largest discrepancy is in the 
statistically precise XcO transition. The lattice systematic errors discussed above may be to 
blame for this, without further calculation we cannot determine this definitively, but we note 
that [?J found that approximately "unquenching" the quark model for this transition did not 
induce a large change in the rate. In Appendix B we consider the effect of using an improved 
vector current and we find that the XcO J/fpl transition is one case where a statistically 
significant change with respect to the local current does occur. This might indicate that this 
transition is particular sensitive to scaling in the lattice spacing. 

The fitted values of j3 in eqn[6]are in rough agreement for all of Xc0,i,2 as one might expect 
for states that in the potential models differ only by small spin-orbit effects. The fitted 
values of A deviate somewhat from the quark-model expected 1 : « : — \ ratio for Xc0,l,2> 
being 1 : 0.62(26) : —0.34(2). Given the approximations and frame-dependence inherent in 
the quark model even such rough agreement is surprising. 

In the Xd,2 J 1^1 transitions, the hierarchy of multipoles expected in the quark model 
(| i±a (0)| > |M2(0)| S> | £3(0) |) is observed. The precise degree of suppression of M2 compared 
with E\ is not easy to compute in non-relativistic quark potential models. However, the single 
quark transition assumption predicts that E^(Q 2 ) = independent of the frame and details 
of the potential and the lattice data do agree with this. Our extracted values of M2(0) or the 
ratio M2(0)/£i(0) all depend upon theoretically undetermined extrapolation in Q 2 and so 
the apparent disagreement with the very small values found in experiment is not yet overly 
concerning. 

The electric dipole couplings for tp' — ► Xco7 and ip" — > Xco7 are, within reasonably large 
statistical errors, in agreement with experimental values and with the "GI" values tabulated 
in [ill 25 . The estimates of Z)-meson loop effects from j3] are comparable to the level of our 
statistical uncertainty. 

For the transition x' c i ~^ J 1^1 we extract a partial width of 21(12) keV which is in 
reasonable agreement with the quark-model estimates of [3], being 14 or 71 keV depending 
upon model details. Since in our study this state has a likely interpretation as the 2 3 P\ 
state of the quark model, without any DD* effects included, this can act as a benchmark for 
models that consider the experimental X(3872) state as being 1 ++ . 

The large £3(0) in the lattice calculation of the x' C 2 ~^ J 1^1 transition is in stark dis- 
agreement with the general result that £3 = for a 3 P2 — * 3 *Si transition. As discussed 
above, this does not appear to be a lattice artifact and so another interpretation must be 
found. Possible explanations include that the x' C 2 state is a hybrid where the gluonic field 
carries spin, but here we concentrate on the most conservative interpretation, that x' C 2 is a 
conventional 3 i ? 2 state. 

An F-wave tensor state ( 3 -F2) is expected to have E\{Q 2 ) = in general from the single 
quark transition assumptions outlined earlier. In addition, a quark model calculation shows 
that M2(0) is highly suppressed: 




25 But note the considerable dependence on the details of the quark-model given in that paper's "NR" results 
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whereas E3 is less suppressed: 



E 3 (\q\ ) oc — ^-exp 



/3 2 ^ V 16^ 2 



with the same leading \q\ 2 dependence as the M2 form factor in a 3 i-2 — ► 3 5i transition (Eqn. 
[TT|) . This pattern that |£/3(0)| > |M2(0)| 3> |i?i(0)| is consistent with the lattice results for 
this x'c2 ~^ J/^Pl transition. These lattice results support the interpretation of this x' C 2 as 
the lightest 3 i ? 2 state having a partial decay width to J/ip'y of 20(13) keV. We are unaware 
of any model calculations of this transition rate. 

The x'c2 ~^ J/i^l transition appears to have reverted to the |£i(0)| > |M2(0)| 3> (0) | 
hierarchy expected for a 3 P2 state and we propose that this state is the 2 3 P2- There is of 
course the possibility that the P and F states are mixed, but this mixing cannot be anything 
like maximal given a large E% appears in only one case and a large E\ only in the other 
case. Within a non-relativistic quark model, the 3 P2 and 3 i*2 states can mix via any tensor 
potential in the Hamiltonian. However in most models this term is suppressed by the inverse 
charm mass squared and so is small. The tensor term, like the hyperfine interaction, is a short 
distance effect and so may not be reliable in these quenched lattice calculations. Because this 
term is suppressed anyway, the errors introduced from this should be small. Another possible 
source of mixing which we lack in the quenched theory comes through DD meson loops - 
such effects are discussed in 0, [32[] ■ 

The partial decay width r(x' c / 2 — ► J/ip'y) extracted from the lattice data is 88(13) keV 



which is comparable with the 53 or 81 keV computed by [18(]. We note that the interpretation 



of excited states in this channel as F and P waves was hinted at in the two-point function 



analysis of 12| where states compatible with being 3 i<2 and 2 3 P2 were found to be nearly 
degenerate. Within the quark potential models used in [l||, the F wave state is expected 
60 — 100 MeV heavier than the P wave. The degeneracy we found may be an artifact of 
"squeezing" these spatially larger states into a 1.2 fm box. 

A discussion of the excited Xc2 states is timely given the observation of a candidate state 



at 3929(5) MeV by Belle in 77 — > DD which has been associated with the 2 3 P2 state 331 ] . 
We lack reliable theoretical estimates of the 77 widths of 3 i*2 and 2 3 P2 states - this calls 



for an extension of the work done in [34| to consider excited states which could use very 
similar operator projection technology to that discussed in this paper. We note here that an 
optimist viewing figure 2 of [iij ] might hope that the statistically insignificant excess at 4080 
MeV could, with increased statistics, become a signal for the other state in the 2 3 i"2 / 1 3 -^2 
pair. 

The J + ( + ) — ► 0~( + -*7 transitions shown in Figure [TBI are consistent with the quark model 
predictions. The electric dipole amplitudes in — > _ ^ + ^7 are much smaller than other 
electric dipole amplitudes since they involve both AL = 1 and a spin-flip. For the magnetic 
quadrupole amplitudes in 2 + ( + ) — > _ ( + '7, the 3 P2 amplitude behaves like Eqn [10] while 
the 3 F2 amplitude is suppressed by a further two powers of |g| which suggests yet again the 
assignment X c2 = 1 3 P 2 , x' c2 « ^F 2 and & « 2 3 P 2 . 



C. Exotic &: crypto-exotic transitions 

The only charge-conjugation allowed transition involving an exotic we compute, rj c \ — > 
J/ip'y, has a rather large partial decay width 115(16) MeV which is dominantly through a 
magnetic dipole transition. Even accounting for the large phase space this is very large on 
the usual scale of magnetic dipole transitions. Conventional cc states can only have magnetic 
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dipole transitions if there is a quark spin-flip and this is suppressed by the large charm-quark 
mass. In a hybrid meson the extra gluonic degree-of-freedom allows an Mi transition without 
spin-flip. This can be seen explicitly within the flux-tube model, where in an Mi transition 
between a conventional L = meson and a hybrid meson the tube absorbs the angular 
momentum 26 . 

The two-point function analysis of [13] suggested that the excited vector state we have 
called ihyb. (having mass around 4.4 GeV) is a crypto-exotic hybrid meson, having con- 
ventional J p = 1 but internally an excited gluonic field. That analysis preferred the 
quarks to be in a spin-singlet in this state. In that case we would expect a large non 
spin-flip Mi transition Y — > r/ c j which is precisely what is seen in our lattice data. Here 
r(V — > 77 c 7) = 42(18) keV which is considerably larger than any other vector to pseudoscalar 
transition. Within the flux-tube model, the non-exotic Y(l^~^ ) and the exotic ?7ci(lhyb) 
(which we find at a mass ~ 4.3 GeV) differ only in being quark spin singlet and triplet re- 
spectively and T(Y — ► n c ~/) « r(?7 c i — > J/ijyy). The lattice data does not strongly disagree 
with this. It would be very interesting to see analogous calculations in other models of excited 
glue to see if this is a general result or one peculiar to the flux-tube model. 

We also consider charge-conjugation violating electric dipole transitions involving exotic 
mesons, 

l-(+) _> 0+(+)7 and 0+M -> 1" (_) 7, where for both we found non-zero signals. 
Roughly speaking, one can gauge the cost of exciting the gluonic field versus exciting conven- 
tional orbital angular momentum by comparing the electric dipole transitions + ( + ) 
having £?i(0) = 0.06(1) and 1~H <-> 0+( +) having #i(0) = 0.127(2). These numbers 
are clearly of the same order, a result that is also true in the flux-tube model 0, Ell - 
whether this is true in other models has not to our knowledge yet been tested. The transition 
+ ( - ) — > 1~( _ ) where -Ei(O) ~ 0.04(1) suggests that L q q = transitions to exotic hybrids are 
also unsuppressed. 

Further evidence for the quark spin triplet nature of the ^cl(l^b) an d the quark spin 
singlet nature of the Y(1^Z ) comes from the fact that Ei(0) for 1~(+) — ► + ^ + ^7 is 0.06(1) 

while for l hj ; b ' — ► + ^ + ^7 it is consistent with zero. The first of these then is S q q = 1 — > 
Sqq = 1 while the second requires a spin-flip, expected to be suppressed by the heavy-quark 
mass. Note that the magnetic dipole transition 1~(+) — > 0~( + ^7 in figure [TT1 is at the scale 
of other "hindered" Ml transitions, consistent with requiring both gluonic excitation and 
quark spin-flip. We note that the spin singlet nature of our Y state does not make it a good 
candidate for the experimental y(4260) state whose potentially large decay rate into ititJ/ip 
would tend to suggest dominance of spin-triplet. 

Here we briefly comment that if this heavy-quark hybrid physics is any guide to the 
behaviour of light-quark hybrid systems, we should expect the GlueX experiment to copiously 
photoproduce hybrid mesons off the meson cloud around a baryonic target. 



V. CONCLUSIONS 

We have demonstrated that the "ideal" excited state operators within a basis of operators 
can be used to successfully extract excited state transitions from three-point correlators. 
Using this technique we have carried out an extensive survey of radiative transitions in 
charmonium with detailed consideration of the phenomenology suggested by the results. 



Details are presented in the unpublished DPhil. thesis of JJD where the following partial widths were 
obtained: r(?7 c i -> J/ipj) w r(Y hyb -> r) c -)) as 30 -> 60 keV 
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Notably we have performed the first lattice QCD calculation of the exotic r\ c \ radiative 
decay and found a large T(r] c i — > J/ip'y) = 115(16) keV. We found statistically significant 
electric dipole and magnetic quadrupole form factors in Xc2 — ► J/ip'Y, calculated for the first 
time in this framework, and have studied excited Xci,2 transitions. Our results for the excited 
Xc2 states suggest that there could be a radially excited 2 3 P2 state and a 1 3 F2 state rather 
close in mass. This signal followed from a clear observation of dominance of E\ over E% 
in one case and £3 over E\ in the other and matches the general expectations of potential 
models. This is relevant given the recent observation of a candidate 2 ++ excited state in 
77 -> DD^. 

Our calculation of magnetic dipole transition widths, such as ip — > 7^7 and hindered 
excited state transitions reflected the expected suppression of the excited state transitions. We 
note that, modulo lattice systematic effects which can be reduced, our method is not troubled 
by the uncontrolled approximations and model dependence inherent in model calculations of 
such suppressed transitions. In light of this we were able to make statements regarding the 
possible influence of closed channel .D-meson loops on the hindered transitions, determining 
that they cannot be as large as suggested in certain studies. We identified a putative non- 
exotic hybrid state having a large magnetic dipole transition amplitude to rj c - the possibility 
that this reflects the non-spin-flip magnetic dipole excitation allowed within the flux-tube 
model (and likely within other models having more than cc content) was raised. 

While the current numerical results may be affected by lattice systematic errors, future 
calculations using the now proven methodology can address these shortcomings by using 
dynamical lattices of sufficient size having a number of different lattice spacings. At some 
point after the introduction of dynamical lattices, the difficulty of dealing with resonant states 
embedded in a multi-meson continuum will have to be addressed. 

We have demonstrated that our results are in agreement with general predictions of quark 
potential models but we are able to go beyond this to make statements about states in which 
there is an excited gluonic field. These results are appropriate for comparison with models 
proposing particular forms for the gluonic excitation. 

The great advantage of this development within lattice QCD as compared to models 
relying upon the non-relativistic motion of quarks, is its immediate applicability to the light 
quark systems. Future efforts will consider photocouplings of light quark mesons, and in 
particular exotics, as these are of central important in the production rates for the GlueX 
project. If the large couplings we find for the 1 ^ state with heavy quarks persist into the 
light quark sector this will confirm the intuition and model results used to motivate the GlueX 
production process. 
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APPENDIX A: DERIVATIVE OPERATORS AT ZERO AND NON-ZERO 

MOMENTUM 

In this appendix we give a few details of the spin-structure we deal with when considering 
two- and three- point correlators. We take the example of the Tl "vector" operators 
considered at p = (000) and p = (100). Using the decompositions in the appendix of we 
find that at zero momentum all the operators we have used have the following behaviour (at 
least in the continuum limit), 

(1— (p = 0,X)\O i \0) oce*(0, A), (Al) 

and there are no non-zero overlaps onto particle of any other J PC (except those lattice 
artifacts that appear suppressed by powers of a). The three-point correlator then involves a 

sum Ea € i (0) A ) £ j(0> A ) = ~9ij- 

At non-zero momentum things are not so simple, with one unit of lattice momentum (e.g. 

p = (100)) we have the following non-zero overlaps in the continuum limit, 

(1— (p, X)\Oi\0) ~ et(p, A), P] e*(p, \) Pi 
<!+-(£ A)|0i|O> ~e ijk e*(p,\) Pk , 
(0 + -{p)\O i \0)^ Pu 

(2 + -(p)|Oi|0> ~ pje^ip, A), Pl e* 00 (p, A). (A2) 

These forms are derived following the decompositions given in Appendix A of [l| 27 . 

So here we have additional contributions from (0, l,2y particles - the (O^)" 1 contri- 
butions we can neglect as these states are very heavy, but the l" 1 we must worry about. 
We can see the entry of such particles in the two-point function spectral analysis; in figure 
[15] we show the spectrum at p = (100) with either i = x, where l" 1 does not contribute, 
or with i = y, z where l -1 can contribute. Also shown are the T x and spectra 
extracted at p = (000) extrapolated to p = (100) using the continuum dispersion relation 
E = ^Jm? + \p\ 2 . We see the entry of the lightest l" 1 state into the T l (p = (100)) spec- 
trum in the case i = y, z but not in the case i = x as expected. We can explicitly exclude 
this state from our three-point analysis by not projecting onto the eigenvector belonging to 
this level. In the results presented in figures 14151 we separately consider the cases i = x and 
i = y, z shown by the triangles and the circles. Note that we did not simply rely upon en- 
ergy matching to determine the levels - a more precise mapping between p = (000) levels and 
p = (100) comes from considering the overlaps, Z. As a concrete example consider the matrix 
elements {l~~(p, A) 1^7^10) = Ze*(p, A) where Z is proportional to the vector decay constant 
and where, for unsmeared quark fields in the continuum limit, Z should be independent of 
momentum p. As can be clearly seen in Table IIV( one can identify the J/tp state with the 
lightest in each case, while the tp' state is the next lightest in two of the cases, but is the third 
state in the case where the l" 1 can contribute. We considered the Z values for the entire set 
of operators when making the state assignments at finite momentum. This approach becomes 
increasingly more cumbersome as the momentum increases so we have not considered it any 
further. 



correcting some minor typographical errors therein 
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FIG. 15: Low lying spectrum extracted from two-point correlators. State assignment (color-coding) 
follows from the consideration of overlap factor as described in the text. 



level 



p = (000) 



p = (100) p = (100) 

Oi = Or, Oi = O. 



0.163(1) 0.173(1) 0.166(1) 

0.190{6) 0.205(10) 0.04(1) 
0.05(25) 0.01(3) 0.201{12) 



TABLE IV: Z values in (1 (p, X)\ipjiip\0) = Ze*(p,X). State proposed to be J/tp in bold, state 
proposed to be ip 1 in italics. 



APPENDIX B: IMPROVEMENT OF THE VECTOR CURRENT 



In [24f | the improvement of the vector current to 0(ma) in a manner compatible with the 
improvement in the anisotropic Clover action was presented. Therein the analysis was for 
heavy-light currents and the renormalization constants were determined perturbatively. Here 
we are interested in heavy-heavy currents and we shall determine Zy non-perturbatively. 
Nevertheless we can consider modifying our local current to include extra terms suggested 
by the improvement scheme to investigate any change in the form-factor values. This can be 
considered to give a crude estimate of how much we might anticipate scaling to the continuum 
to affect our results. 

Our Clover action uses r = 1 and as such our local vector current is not automat- 
ically improved - at tree level the improved current is given by where ^/(x) oc 

^1 + a s di7j D j \ i/j(x) is the "rotated" field used in the construction of the improved action. 

The improvement parameter d\ = ^(1 — £r) + 0(rriQat) ~ —0.5 for renormalised anisotropy 
£ = 3.0. The quark mass parameter that appears in our Clover action, moat has the value 
0.0401 which is clearly small, while the mass in spatial lattice units is three times as large, 
but still might be argued to be small. 

We have attempted to construct the improved current by using the tree-level equation of 
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motion 

(a t 74 V 4 + a s ^jjVj + m a t ^ ip = Q (Bl) 

to eliminate the derivative in the expansion to 0(a s ) of the improved current ^ r y fJi ^. This 
yields the result 

^74^ oc Vh^V* ~~ dia 8 dj (ipajA'tf)) (B2) 

# 7i * oc (1 - 2m a i |)^ 7 ^ - d^d^a^) (B3) 

where a common proportionality factor is ignored given that we determine Zy non- 
perturbatively using the pseudoscalar form-factor at Q 2 = 0. At q = (000) the second term 

does not contribute and we would expect Zv^ t) = ^~ a * 77^1 ~ 1-H for the parameters in 
our action. Note that in section [III Al we found 1.11(1) for this ratio using a non-perturbative 
extraction. 

We computed form-factors using the improved currents for a few of the transitions con- 
sidered in this paper. The only effect of a considerable size was found in the XcO ~~ > J/^l 
transition as shown in figure [ToTa) . Note that the addition of the improvement brings the 
Clover data into better agreement with the DWF data on the same lattices as we might ex- 
pect given the automatic 0(a) improvement one has with DWF. We note that for the same 
correlators projected on to the excited state ij/ such a large difference with respect to the local 
current was not seen, see figure [ToTb). We were unable to find any other large effects due to 
improvement, e.g. consider J/ip — > i] c -f shown in figure [T71 These observations may suggest 
that except in certain notable cases (the scalar), the discretization errors on our results are 
relatively small. 
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FIG. 16: Scalar- Vector E\ transition form-factors using local and O(ma) improved current with 
anisotropic Clover quarks and the fit from a study using the local current with DWF fermions on the 
same anisotropic lattices, (a) \c0 ~ * J/ipl, (b) ip' ~ * Xco7- 
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FIG. 17: J/ij) — > r/cj transition form-factor using local and 0(ma) improved current with anisotropic 
Clover quarks and the fit from a study using the local current with DWF fermions on the same 
anisotropic lattices. 
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